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ABSTRACT 

A three dimensional model for the annual and semiannual variations of the 
thermosphere is presented in which energy and diffusive mass transport associ- 
ated with the global circulation are considered in a self consistent form. It is 
shown that these processes play a major role in the thermosphere dynamics and 
thus account for a number of temperature and composition phenomena such as 

(a) The helium and oxygen bulges in the winter hemispheres at lower alti- 
tudes, 

(b) The relatively large temperature variations (~ 26%) in the annual com- 
ponent exceeding by a factor of three those inferred from the satellite 
drag data and which would be much too large for a thermosphere in dif- 
fusive equilibrium, 

(c) The 7%-variation of the exospheric temperature at the poles in the 
semiannual component associated there with a significant depletion of 
oxygen and helium during the equinox, 

(d) The nearly height independent and relatively small amplitude in the 
total mass density which is responsible for the fact that the latitude 
dependence of the semiannual effect has so far eluded observations 
from the satellite drag data. 
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THEORETICAL MODEL FOR THE LATITUDE DEPENDENCE OF THE 
THERMOSPHERIC ANNUAL AND SEMIANNUAL VARIATIONS 


1. INTRODUCTION 

The seasonal variations in the thermospheric density, as deduced from 
sateUite drag data (Jacchia et al., 1966), have been of great concern for a long 
time. It should be expected that the annual and semiannual components are in a 
relation that corresponds approximately to the corresponding harmonic compo- 
nents of the solar heat input. Thus, energetically, a major problem emerged 
when it was observed that the semiannual effect appeared to be latitude inde- 
pendent (Cook, 1970) and that its magnitude seemed to be equal to that of the 
annual component. 

The annual component of the thermosphere has so far been observed pri- 
marily from satellite drag data which were interpreted in terms of temperature 
variations (e.g. Jacchia and Slowey, 1967; Volland and Mayr, 1971). However, 
the existence of the winter anomaly in the F 2 region has led very early to the 
speculation (Johnson, 1964) that the annual component in the thermospheric wind 
field would decrease the oxygen content in the summer hemisphere, a hypothesis 
that was confirmed by Mayr and Volland (1971). The discovery of the winter 
helium bulge (Keating and Prior, 1968) brought further evidence for the significance 
of the large scale circulation to which it was attributed by Johnson and Gottlieb 
(1970). 

Until recently all models have attributed the semiannual density variations 
to corresponding variations in the thermospheric temperature (Jacchia et al., 
1966; Volland, 1969; Johnson, 1964; NeWell, 1968). 

Jacchia et al., (1966) associated temperatures to the observed density vari- 
ations, Volland (1969), in following a suggestion by NeweU (1968), explained the 
semiannual temperature amplitudes by waves leaking from the mesosphere up 
into the thermosphere and Johnson (1964) proposed annually varying large scale 
meridional circulations which subtract more energy from the thermosphere 
during solstice than during equinox. 

Cook (1967) proposed changes in the turbopause level to explain the semi- 
annual variations at higher altitudes (800 km) where helium becomes the major 
constituent. This mechanism, however, which invokes variations in the eddy 
diffusion coefficient below 120 km, and consequently would also change the rela- 
tive concentration of atomic oxygen was examined by Harris and Priester (1969) 
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and Cook (1969) who concluded that it would be inconsistent with the semiannual 
variations in the thermosphere up to 700 km. In analyzing the semiannual effect 
in the F 2 region Mayr and Mahajan (1971) concluded that there should be signifi- 
cant semiannual variations in the (0)/(0 2 ) and (0)/(N 2 ) ratios in the lower 
thermosphere with maxima occurring during equinox an "observation” that ap- 
peared to be consistent with the few reported rocket measurements. In explain- 
ing these variations in the neutral composition a circulation model was proposed 
by Mayr and Volland (1971). In this model the semiannual component of the 
meridional wind circulation was attributed to a primary energy source at high 
latitudes (presumed to be associated with the semiannual component in the oc- 
currence of magnetic stroms) and to a much weaker secondary source at low 
latitudes (associated with the solar EUV heat input). Under the influence of this 
circulation, which was matched to agree with mesospheric wind measurements, 
atomic oxygen is redistributed semiannually such that during equinox O is de- 
pleted at high latitudes (and to a much lesser degree at the equator) and enhanced 
at mid to low latitudes. Since this mechanism produces during equinox and at 
mid to low latitudes enhancements in the concentration of atomic oxygen without 
'invoking temperature variations, it was concluded that the models that relate (at 
these latitudes) the semiannual density variations solely to temperature effects 
would probably overestimate the semiannual temperature amplitude. In a recent 
paper, Jacchia (1971) completely abandoned the temperature derivation from the 
semiannual effect in his model, stating that all difficulties in the description of 
the semiannual density variations could be removed if they were not related to 
temperature variations. Furthermore, Brinton and Mayr (1971) concluded from 
the small semiannual component in atomic hydrogen at mid latitudes that the 
semiannual variations in the exospheric temperature should be negligible there. 

In comb inin g all these evidences one might therefore suggest that a signifi- 
cant part of the semiannual density variations is associated with changes in the 
neutral composition. 

Semiannual variations in the neutral composition, produced by global circu- 
lation must however be accompanied by a semiannual, latitude dependent tem- 
perature component. 

It will be the subject of this paper to discuss the latitude dependent aspects 
of the annual and semiannual variations quantitatively on the basis of a theoretical 
model. 


2. THEORY 

Figure 1 shows a very simplified block diagram for the thermosphere dy- 
namics. It illustrates the links between the composition and gas temperature 
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Figure 1. Schematic diagram for the physical processes that influence the temperature 

and composition of the thermosphere. 

on one hand and the various physical processes and atmospheric parameters on 
the other hand. 

Suppose the heat input Q is known as a function of time, height, and latitude. 
Some of this energy input is conducted down into the lower atmosphere, thus af- 
fecting the temperature distribution and with that the composition. 

However, a second and not less important energy train goes through the 
global circulation which is significantly affected by ion collisions and thus by 
the ion density. In this mechanism, which is explicitly neglected in one dimen- 
sional models, energy is exchanged between summer and winter through adiabatic 
expansion and contraction and this in turn affects the temperature and the 
thermospheric composition. 

In parallel, the thermospheric circulation induces diffusion which depends 
upon the eddy coefficient, and this diffusion process affects the composition 
directly. Any variations in the composition are associated with variations in the 
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total pressure field and this in turn influences the global circulation and with 
that adiabatic heating and with that again temperature and composition. So it is 
apparent that all the processes that have been surrounded by the dashed box in 
Figure 1 are interconnected. We shall show in this paper that the link between 
diffusion and energetics is extremely important for seasonal variations in the 
thermosphere. 

To incorporate the diffusion processes in a self-consistent form, a two 
component model is adopted for the three dimensional thermosphere dynamics. 
With the subscripts 1 and 2 referring to the two atmospheric constituents, the 
equations of conservation of mass, momentum and energy are 


^ P\ 2 - 

— " °i + L i + 3-7^1 w i) + div s (P 1 Vj) - 0 


(la) 


^ Pi 2 

- ° 2 + L 2 + __ ( p 2 w 2> + div 0 (p 2 V 2 ) = 0 


(lb) 


P. 


B P1 1 

f D nij + K m 

Br +l 

V m ! (D + K) 

3 p 2 , 

m 2 + K m 

B r ( 

. m, (D + K) 

\ t 

B Vl 
B t 

'l V j + 0 ^ 2 

bv 2 

T* ± 


gP2+ 4> Pi p 2 ( w 2 - Wl ) = o 


(2a) 


(2b) 


3 2 V 

2 1 1T + v 2 v 2 + 4 Pi <v 2 - + 2 n cos fiu 2 f-i) — - = 0 ( 3b ) 

B r 2 


BU. 


P. 


. * • . B 2 U. 1 3p, 

bT +Vl sin 1 u i + $p 2 < u i - u 2 ) - 2 n cos ev^-v — - +- jj = ° ( 4a > 

^ J B r 2 r 


BU„ 


B 2 U„ 1 B 


p2 1tF + I ' 2 sin * IU 2 + <Pp t (U 2 - Uj) - 2 n cos 9 V 2 j> - 77 2 +1 —5 = 0 ( 4b ) 


B r 2 r 
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2 X^rrij m 2 J 3 t 
with 



3Wj 3 W 2 

— + P 2 — + Pi div 0 u i + P 2 div 0 U 2 =Q(5) 


and 


div g 


1 

sin 6 


A. (sin 9) 


( 6 ) 


p v p = mass densities 

Wj, W 2 = vertical transport velocities 

Vj, V 2 = horizontal transport velocities in the longitudinal direction 

Uj, U 2 = horizontal transport velocities in the latitudinal direction 

r, 6 - radial distance, colatitude 

°i’ °2 - chemical production terms 

Lj, L 2 = chemical loss terms 

V v v 7 ~ collision frequencies with ions 

<{> = k T 

(D + K) m x p x 

k = Boltzmann constant 
T = temperature 

D, K = molecular eddy diffusion coefficient 
p r p 2 = partial pressures 
r) = viscosity coefficient 
k - conductivity coefficient 
Q = heat input rate 
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I = dip angle 

n = angular frequency of the earth's rotation. 

In Equations (1) we have assumed the chemical production and loss terms to 
be zero in this paper, since we shall emphasize the dynamic aspects. In a sub- 
sequent paper we shall discuss the influences of 0 2 dissociation, recombination 
and the transport of chemical energy. 

Equations (1) through (5) are solved under the following assumptions: 

(1) The variables are separated into time t, latitude and altitude r com- 
ponents in the form 


_ j co( t - t „) 

P = P 0 + p (O Pp (*) e P 1 

j o) ( t “ t 

P = P 0 + P (r) P p ( 9 ) e 
T = T 0 +T(r) P T (6) e J “ ( 4 “ t T ) 

(7) 

W = W ( r) P ff ( 9 ) e J “ (t " tw) 

V = V (r) P y ( 0 ) e , “ (t_t v ) 

U = U(r)P u (0) e J “ (t_tu) 


where oj are the angular frequencies corresponding to the periods of one and 
one half years, and tj are the phases (times of maxima) for the individual 
atmospheric parameters. For the latitude dependences, P , P p , etc., expansions 
into spherical harmonics are adopted thereby neglecting higher order terms: 


p — P — P — P — 
P P 


Pj ( 6 ) : annual 


P ( 6 ): semiannual 


(8) 
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with 


sin 0: annual 

sin#: semiannual 

sin 6 cos 8 ; annual 


P o = 


sin (9 cos 0: 


semi annual 


Pj = cos 8 


P 2 = 1/2 (3 cos 2 6 - 1) 


(9) 

( 10 ) 


( 11 ) 


being Legendre polynomials. This simplification is made for mathematical con- 
venience and in order to keep the computation time in limits that guarantee the 
usefulness of the model as a tool. It can be shown however that these functions 
approximate rather well the corresponding fundamental planetary waves at 
thermospheric heights (Volland and Mayr, 1971). 

(2) Considering the functional forms (6) through (11) in the continuity equa- 
tions (1) through (5), perturbation theory is applied. This leads (after differen- 
tiating some equations) in a straight forward manner to a set of 9 complex or 18 
real second order differential equations which will not be shown here in the inter- 
est of conciseness. These equations describe the height dependence in the am- 
plitude and phases of the atmospheric parameters. 

The numerical scheme in performing the integration is essentially that of 
Lindzen and Kuo (1969). It will be seen in section 4 that in the case of the annual 
wave the amplitudes of N and of He become quite large so that the validity of 
perturbation theory may &e questionable. However these amplitudes become 
greater than one primarily in those heights where N 2 and He are the minor con- 
stituents. Therefore they cannot influence significantly the energetics and dy- 
namics of the total gas. Amplitudes greater than one are of course not realistic 
in a linear theory and should be considered as upper limits. 

The boundary conditions employed in the calculations where the foilwing: 

(1) At the lower boundary it was assumed that the characteristic time of the 
atmosphere is so long that no significant variations can be excited with the 
seasonally varying energy input. Accordingly we assumed 
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p = p = T = W = U = V= 0 at 80 km. 


(12) 


This assumption is strictly speaking somewhat arbitrary, and therefore the 
solution is meaningless at the lower boundary level. It can be shown however 
that at higher altitudes above 120 km the solution can be considered as unique. 

(2) Because of the low internal energy and heat input rate it was assumed 
that at the upper boundary the transport of mass, energy and momentum have 
negligible effects upon the thermosphere dynamics. It was therefore postulated 
that 


|-^ = pL= JI = ° at 500 km. (13) 

or or 3 r 


Again it must be noted that this assumption is actually not completely valid. The 
effect from this ambiguity, however, becomes negligible below 400 km. 


3. INPUT PARAMETERS 

For the heat input due to the solar radiation a height dependence 

* / f — j- ^ /w 

Q =? Q 0 e 0 ( r ^ r Q - 1 20 km) (14) 


is adopted. Choosing a scale height of H = 50 km we assume the same distribu- 
tions that has been used in our model for the diurnal variations in the neutral 
composition (Mayr and Volland, 1972). 

From a harmonic analysis of the solar heat input into the atmosphere 
(Volland and Mayr 1971) the amplitudes of the annual and semiannual components 
are approximately related to the diurnal component Q° through the relations 


Q 0 a = 0.4 
Qf = - 0.075 Q° 


(15) 
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where Q A and Q® A are the coefficients for the latitude dependences P°(0) and 
P° (6) in the annual and semiannual heat input distributions, respectively. 

The actual heat input distribution in the thermosphere is quite uncertain. 
For this reason Q ° has been chosen as a free parameter which was determined 
in Mayr and Volland (1972) to 


Q° = 2 x 10 -7 ergs/cm 3 sec at 120 km 

such that the resulting temperature amplitude agreed with observations. Cor- 
respondingly we shall choose in this paper the values 

Q a = 8 x 10 8 erg/ cm 3 sec (16a) 

Q 8A = - 1.5 x 10 _8 erg/cm 3 sec. (16b) 


The input (16a) will be the only source for the annual variations in our model. 
For the semiannual component, however, the source (16b), which is due to the 
sun's semiannual migration across the equator, is only one contribution. 

In addition to (16b) at least the contribution from the semiannual component 
in the occurrence of magnetic storms must be considered. 

For the electron density, which enters into the collision frequencies v and 
v 2 an altitude distribution of the form 


N = N e 


<r- r ) /H( r) 


(17) 


is used, with H(r) varying in a form 


H (r) = hj (r - r m ) + h 2 (r - rj* 


such that the scale height becomes 
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Hj = 80 km at 120 km. 


and 


H 2 = - 150 km at 400 km. 


For the F 2 maximum electron density a value of 


N = 6 x 10 5 /cm 3 

m ' 


(18) 


is adopted at r m = 300 km. With the neutral-ion drag coefficient of 


€ - 7.5 x 10 10 cm 3 / sec (Dalgarno, 1964), 


the collision frequency = v 2 = v is derived from 


" = N e 



(19) 


where a> L = 1.15 x 10 2 /sec is the Larmor frequency for ions and n is the time 
average number density of the neutral atmosphere. The electron density is as- 
sumed to be latitude independent, and the expressions sin 2 I U x 2 in (4a) and 
(4b) are approximated by 


sin 2 I Uj 2 ~ 


0.413 U 


1,2 


1 ^ 0.69 


U 


1 , 2 - 


annual 


semi annual 


( 20 ) 


For the Coriolis force in (4a) and (4b) the following approximations were neces- 
sary in order to facilitate the separation into spherical harmonics: 


cos 2 e ~ 0.333 P°: annual 

(21) 

cos#pO~0.4 P j : semiannual 
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The heat conductivity is computed from the function 


k = 5.25 x 10 3 T 1/2 /M 


( 22 ) 


with M being the mean molecular mass in AMU. This form is in close agree- 
ment with the conductivity coefficients for the individual atmospheric constituents 
quoted in Harris and Priester (1962). 

For the viscosity coefficients a form 


rj = t ?0 T" 2 (23) 

is adopted with a value of rj Q = 1.3 x 10" 5 from Nicolet (1960). 

Because relatively little is known about the eddy diffusion, a height independent 
average value of 


K = 4 x 10 6 cm 2 /sec 


is taken for the diffusion coefficient from Colgrove et al., (1965). However, this 
value will be modified in order to illuminate the sensitivity of the solutions to 
the eddy diffusion. 

For the molecular diffusion coefficient the function 


D = D 0 (T/T s ) 1,7S (p/p s ) _1 (24) 

is used where T s and p s are standard temperature and pressure, p being the 
sum of the partial pressures for the major constituents. 

From the work of Walker (1961) a value of D Q = 0.26 was taken for the dif- 
fusion of oxygen through 0 2 and N 2 . Considering that the diffusion coefficient 
depends on the square root of the reduced mass for the collision partners, a 
value of D he = 0.50 is adopted for helium diffusing through O and N . 
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For the time average temperature distribution above 120 km a Jacchia (1964) 
model is adopted with an exospheric temperature of = 1050°K, while below 120 
km the temperatures are taken from CIRA (1965). For the densities at 120 km 
the values are taken from Jacchia ( 1964 ). 

It is assumed that the time average densities are in diffusive equilibrium. 
Since over the entire altitude range a constant eddy diffusion coefficient is used 
the transition between mixing and diffusive separation is continuous in our model. 


4. ANNUAL VARIATIONS 

Figure 2 shows in schematic form the basic elements of our theoretical 
model for the annual variations. The solar heat input maximum in the summer 
hemisphere causes the gas to expand at the summer pole. As a result, a circu- 
lation cell forms in which kinetic energy is transported from the summer to the 
winter hemisphere thus damping the temperature amplitude that would develop 
in the absence of this advective energy loss. At the same time this thermospheric 



Figure 2. Schematic picture for the latitude dependence of the annual component of the various 
physical parameters with the important processes emphasized in circles. 
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circulation carries the minor and lighter constituents O and He from the summer 
to the winter hemisphere thus causing a redistribution in the thermospheric 
composition which is characterized by the He and O bulges in the winter hemi- 
sphere. 

Figure 3 shows the amplitudes for N 2 , O, He, and T g at the summer pole 
from a solution of our theoretical models. With Am / 0, our N 2 - O model is 
identified. In this model the diffusion of both these constituents is considered 
and thus it represents also a realistic condition for the energetics of the 
thermosphere which is essentially determined in the region where N 2 and O 
are the major constituents. With Am = 0, our He- model is characterized in 
which the diffusion of He through N 2 and O is treated, assuming that the mean 
molecular mass of these constituents combined does not vary as a function of 
time. 


ANNUAL VARIATIONS 



Figure 3. Relative amplitudes for the 0-, He-, and Hj-concentrations, the mass density p and 
the gas temperature Tg in the annual component. The values are valid for the summer pole and 
are the coefficients of the spherical function P^(0) = cos 6 {0 = polar distance) which describes 
the latitude dependence. Positive signs indicate positive, maximum values, negative signs in- 
dicate negative, minimum values in the northern summer hemisphere. Solid lines are taken from 
the Nj-O-model which describes correctly the energetics and thus the gas temperature. Dashed 
lines stem from the He-model which is energetically unrealistic corresponding to the assump- 
tion of diffusive equilibrium for the major constituents Nj and 0 up to 500 km. Note the large 
difference between the “diffusive equilibrium” solution and the “diffusion” solution for the gas 
temperature Tg as well as the small and nearly height independent amplitude of the mass den- 
sity p. 
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With plus signs, it is indicated that a particular atmospheric parameter is in 
phase with the solar heat input or that it increases in the summer and decreases 
in the winter hemisphere, while the minus signs indicate a decrease in the sum- 
mer and an increase in the winter hemisphere. 

From this solution it is then apparent that a large winter helium bulge 
develops throughout the entire altitude range above 100 km. This result is 
basically in agreement with (a) the observations of Keating and Pior (1968) who 
detected the winter helium bulge in the satellite drag data and (b) the theoretical 
results of Johnson and Gottlieb (1970) who adopted a more or less reasonable 
wind field for their calculations. 

The relative temperature amplitude in our He solution is about 4% and thus 
very close to the 5% variation inferred from the annual component in the satel- 
lite drag data (Jacchia, 1965). This kind of agreement is very si gnifi cant be- 
cause our He solution corresponds essentially to one of diffusive equilibrium for 
the major constituents up to 500 km, a condition also presumed to be valid in 
Jacchia's model. However, deviating from Jacchia's model which assumes that 
the temperature variations are zero at 120 km, our results show there a rela- 
tively large temperature amplitude. And this result had to be expected consider- 
ing that the response time of the thermosphere is very much shorter than the 
period of a year even within the lower thermosphere, a situation quite different 
from those encountered in the diurnal and magnetic storm variations. 

With solid lines the solution from our N 2 - O -model is represented. Since 
in this model the diffusion of the major constituents O and N 2 is considered in 
a self consistent form, it provides a realistic picture for both the O- and N 2 ~ 
composition as well as for the gas temperature. 

This solution shows the formation of a winter oxygen bulge which is par- 
ticularly pronounced within the lower thermosphere. It confirms thus the results 
by Mayr and Volland (1971) which were based on a wind field inferred from meteor 
trail observations. With a factor of three increase in O between summer and 
winter pole at 120 km, the formation of the winter anomaly in the F region can 
thus easily be tied to the global circulation as has been suggested by Johnson 
(1964) and King (1962). 

Although the diffusion velocities are generated within the entire thermo- 
sphere, they are effective primarily within the lower thermosphere and meso- 
sphere where the ambient density is relatively high. Thus at higher altitudes 
the thermal expansion, associated with the temperature peak in the summer 
hemisphere, becomes increasingly more important. The consequence is that in 
the summer hemisphere the depletion of atomic oxygen, associated with the global 
circulation, becomes more and more compensated by the thermal expansion at 
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higher altitudes. This is apparent in Figure 3 where the oxygen distribution 
exhibits an isopicnic level at exospheric heights. Below this level, the diffusion 
process dominates, and an oxygen bulge develops in the winter hemisphere; 
above it, the oxygen density is in phase with the temperature and thus peaks in 
the summer hemisphere. In helium, which is even more sensitive to the diffu- 
sion process and relatively insensitive to the temperature variations, such an 
isopicnic level can not develop and thus the winter helium bulge remains intact 
throughout the thermosphere. 

One of the major results that comes out of our model is that the diffusion 
process has a very strong effect upon the temperature amplitude. With the de- 
pletion and enrichment of oxygen in the summer and winter hemispheres, re- 
spectively, (induced by the diffusive redistribution) the total pressure variation 
in the annual component is substantially damped. Hence also the meridional 
velocities are significantly reduced and with that the adiabatic energy exchange, 
the consequence being that a much larger temperature variation can be main- 
tained. In a way the diffusion process constitutes an effective viscosity which 
damps the meridional velocities and thus tends to isolate the solar heat input. 

In the N 2 ~0- model (solid line) which is appropriate for describing the tempera- 
ture distribution, the amplitude is about 13% at exospheric heights thus about 
three times larger than that from the He-model (dashed line) in which the dif- 
fusion of the major constituents has not been considered. This temperature 
amplitude is also by about a factor of three greater than the temperature ampli- 
tude in the annual component of Jacchia's model. 

A temperature amplitude of this magnitude, if associated with a diffusive 
equilibrium model, would produce meridional velocities of such magnitude that 
the corresponding energy exchange would essentially wipe out the temperature 
variations. It must be emphasized that such a large temperature amplitude can 
only be understood in the framework of a model in which diffusion and adiabatic 
heating are considered in a self-consistent form. 

Variations in the total mass density are finally shown as a thick solid line 
in Figure 3, the plus sign indicating that, as expected, the density increases in 
the summer hemisphere. This quantity which is essentially observed in the 
satellite drag data reflects the superposition of two effects that tend to cancel 
each other: the thermal expansion of N 2 and the depletion of O in the summer 
hemisphere (and vice versa in the winter hemisphere). In this respect the mass 
density has thus a characteristic quite similar to that of the total pressure vari- 
ation which determines the wind field. It is apparent then from Figure 3 that the 
amplitude in the density is relatively small with some tendency even to decrease 
between 260 and 350 km, the latter reflecting upon the increasing relative abun- 
dance of atomic oxygen at these heights. 
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Altitude (km) 


Although the details in the height distribution of the mass density depend on 
a number of factors to be discussed in the following, it is characteristic for its 
amplitude that it tends to remain constant with height. An interpretation of the 
satellite drag data in terms of a diffusive equilibrium model would therefore 
lead to the conclusion that this height structure is produced by a wave that travels 
from the lower atmosphere up through the thermosphere; while in reality we just 
see the manifestation of a highly variable composition which in fact is not in 
diffusive equilibrium. 

In Figure 4 the annual component of the thermospheric wind field is presented 
from the N 2 -Omodeb The maximum vertical velocities, W, are directed upward 
and are shown for the summer pole, the maximum meridional velocities, U, blow- 
ing toward the winter hemisphere are shown at the equator, and in the northern 
hemisphere the zonal components, V, are eastward in winter and westward in 
summer with their maximum magnitudes taken at 45° latitude. The zonal wind 



W Icm/sffl U,V (m/sec) 


Figure 4. Amplitudes from the thermospheric wind field in the annual component of the Nj-O- 
model. The maximum vertical velocities, W, are directed upward and are shown for the summer 
pole, the maximum meridional velocities, U, blowing toward the winter hemisphere are shown at 
the equator, and in the northern hemisphere the zonal winds, V, are eastward in winter and west- 
ward in summer with their maximum amplitudes at 45° latitude. Dashed lines are for the O-com- 
ponent, solid lines are for the (^-component. Note that the difference between both constituents 
is particularly large for W, but is also not negligible for U. 
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which is due to the geostrophic component of the wind field decreases toward 
higher altitudes thus indicating the diminishing influence of the Coriolis force. 

The meridio nal wind velocity increases with height reflecting the increasing im- 
portance of ion drag and viscosity which produce winds parallel to the pressure 
gradient. The magnitudes of both wind components are consistent with the re- 
sults from a three dimensional theoretical model (Volland and Mayr, 1971). 

It is apparent that the horizontal transport velocities for O and N 2 can be 
different by as much 25 % at exospheric heights. Therefore the diffusion effect 
of the horizontal wind field can not be entirely neglected as this is possible for 
the diurnal variation. 

Much more significant is the difference between the vertical transport velo- 
cities of both constituents with the oxygen velocity being by about a factor of two 
larger than that of N 2 . It is this difference which accounts for the decrease in 
the scale height of atomic oxygen during summer and its increase during winter. 

The eddy diffusion properties of the mesosphere are not well known. For 
this reason it is of interest to discuss the influence of eddy diffusion upon the 
thermospheric structure. In Figure 5 two solutions are presented from the N 2 ~ 
O- model for eddy diffusion coefficients of 2 x 10 5 6 and 8 x 10 6 cm 2 /sec. From 
Figure 5 it is apparent that all atmospheric quantities are significantly affected. 
The most drastic effect occurs in the amplitude of the oxygen concentration which 
increases by almost a factor of two at 120 km height when the eddy diffusion co- 
efficient has increased by a factor of four. Since the relative variations of O are 
primarily responsible for the temperature amplitude as pointed out before, the 
increase of the eddy diffusion coefficient also increases the temperature ampli- 
tude significantly. This in turn leads to a stronger compensation of the oxygen 
depletion by the superimposed thermal expension with the final effect that at 
higher altitudes the variation of the total density due to variations in the eddy 
diffusion becomes damped. 

Similar effects are apparent when the total content of atomic oxygen is 
varied. E.g., an increase of (O) at 120 km from 7.5 x 10 10 cm -3 to 1.5 x 10 11 
cm -3 leads to a substantial reduction (-30%) of the relative variation of O and 
T . With the uncertainties in the eddy diffusion coefficients and in the oxygen 
concentration the model is thus far from being unique and is therefore only 
considered a means of identifying some physical aspects of the thermosphere 
dynamics. 


5. SEMIANNUAL VARIATION 

A number of sources have been proposed for the excitation of the semiannual 

effect in the thermosphere. 
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ANNUAL VARIATIONS 



RELATIVE DENSITY VARIATIONS AN/N 0 RELATIVE TEMPERATURE VARIATIONS AT/T 0 


Figure 5. The Nj-O-model shown for different eddy diffusion coefficients. 

Note that all parameters are affected including Tg. 

(a) the non-linear energy coupling between the annual and semiannual 
variations , 

(b) the solar EUV-component (Eq» 15) which is in a certain relation to 
the annual heat input, 

(c) the auroral heat input associated with the semiannual component in 
the occurrance of magnetic storms and 

(d) the energy coupling from the lower atmosphere where the semiannual 
effect is relatively strong. 

We performed calculations on the non-linear coupling from the annual to the 
semiannual variations associated with adiabatic cooling and found it to be neg- 
ligibly small at least for the latitudes dependent part of the semiannual variation. 
This result may be surprising in view of the relatively large variations in all 
atmospheric parameters of the annual component. However in the energy coupling 
the sum of both density components O and N 2 enters, and as shown in section % 
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these quantities tend to cancel each other as is apparent from the small ampli- 
tude in the mass density in Figure 4. 

It has been suggested by Mayr and Volland (1971) that a latitude dependence 
in the semiannual variation of the temperature and composition should be pri- 
marily excited by the magnetic storm component at high latitudes and by a much 
weaker source at low latitudes associated with the solar EUV heat input. Recent 
observations by Hedin et al., (1972) confirm this concept except with regard to 
the low latitude structure which could very well be masked by the dominance of 
the height latitude source. 

We shall adopt here the simplified model that the latitude dependence in the 
semiannual variation is in a first approximation excited by a source of the form 


Q sa = Q 0 sa P 2 ( 9 ) 


which peaks at the poles during equinox reflecting the dominance of auroral 
heating. Since the magnitude of the magnetic storm component is not known, 

Q 3a is considered as a free parameter to be determined by matching the theory 
with the exospheric temperature derived from the data of Hedin et al., (1972). 
This value is 


Q aa = 2 x 10~ 8 erg/cm 3 sec at 120 km height, 


which is of the same order as the heat input due to EUV (Eq. 16b), however re- 
versed in sign. For the height dependence of Q sa we adopted the same formula 
as for in (14) with H = 50 km. 

The results of our model calculations are shown in Figure 6. Again with 
dashed lines the He-model and with solid lines the N 2 -0- model is represented. 
The thick solid line gives the total mass density of the N 2 -O model. The ampli- 
tudes are the coefficients of the P 2 -function. They therefore describe the actual 
values at the poles. According to the dependence of P 2 on colatitude 6 (see Eq. 
11) plus-signs indicate that a particular quantity is enhanced at the poles and 
decreased at the equator during equinox. Minus signs indicate a decrease at the 
poles and an increase at the equator. 

The features apparent from Figure 6 are completely analogous to those for 
the annual component in Figure 3. The thermospheric circulation that is pro- 
duced by the equinoctial pressure and temperature maxima at the poles induces 
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SEMIANNUAL VARIATION 



Figure 6. Relative amplitudes for the physical parameters of density and temperature of the 
semiannual variation. The values are valid for the poles andarethe coefficients of the spherical 
harmonic P 2 (0)- Positive signs indicate at the poles maximum and minimum values during 
equinox and solstice, respectively. Negative signs indicate the opposite. The N 2 -0-model and 
the He-model are differentiated with solid and dashed lines, respectively. Note the large dif- 
ference between “diffusive equilibrium" solution (dashed line) and “diffusion" solution (solid 
line) for Tg, and the height independent and very small amplitude of the mass density p. 

there a depletion in the population of oxygen and helium. This effect is pri mar ily 
pronounced at low altitudes (minus signs for oxygen and helium). It is graduately 
compensated by the thermal expension in the distribution of O such that an isopic- 
nic level develops near 450 km with oxygen being in phase with the temperature 
at. higher altitudes. 

In helium the thermal expension is less effective, and an isopicnic level 
cannot develop up to very high altitudes due to the more pronounced depletion 
effect. Therefore He decreases throughout the termosphere in the region of 
enhanced temperature. All these features are observed by Hedin et al., (1972). 

The difference between the temperature amplitudes from the He-model 
(dashed line) and the N 2 -0- model (solid line) is again very large and thus re- 
flects upon the importance of the diffusion process for the thermospheric 
energetics in the seasonal variations, an effect that has been discussed in detail 
for the annual component. Since N 2 and O are the major constituents up to 500 
km, the temperature amplitude from the N 2 -O- model (solid line) is considered 


20 



as realistic. It is this quantity which has to be compared with the temperature 
derived by Hedin et al., (1972). The distribution of the total mass density is 
flat and its amplitude is very small similar to the annual variations. The latter 
accounts for the fact that the satellite drag data have so far not revealed a lati- 
tude dependence in the semiannual variation as predicted by Mayr and Volland 
(1971). 


6. CONCLUSION 

A three dimensional model of the thermosphere is presented in which energy 
and diffusive mass transport associated with the global circulation are considered 
in a self consistent form. It is shown that these processes play a major role in 
the thermosphere dynamics and thus account for a number of composition and 
temperature phenomena which are now being observed. In particular, the model 
describes 

* 

(1) In the annual component 


(a) the helium and oxygen bulges in the winter hemisphere, 

(b) the relatively large temperature variations (26% between winter 
and summer poles) which exceed by a factor of three the temper- 
ature variations from the Jacchia model and which would be en- 
tirely impossible in a diffusive equilibrium model, 

(c) the nearly height independent amplitude of the total mass density 
which is in our model the manifestation of a highly variable com- 
position. 

(2) In the semiannual component 


(a) the significant latitude dependence in the composition with a de- 
pletion of oxygen and helium at the poles and a corresponding 
enhancement in the density of these constituents at the equator dur- 
ing equinox, 

(b) the latitude dependence of the temperature with a 7% variation at 
the poles where the temperature peaks during the equinox, 

(c) the negligibly small amplitude in the latitude dependence of the 
mass density which is responsible for the fact that the latitude 
dependence of this effect has so far eluded observations from 
satellite drag data. 
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It is the self consistent treatment of the diffusion and energy transport 
mechanisms that accounts for the agreement between our model and the latitude 
dependent components of the seasonal variations now being observed. Although 
sophisticated with regard to the physics of the thermosphere dynamics, the 
model is mathematically however relatively primitive. The application of per- 
turbation theory and the restriction to the first harmonics with their simplifica- 
tion in the description of ion drag and Coriolis force characterize the approxi- 
mations. This becomes particularly evident in our calculations when amplitudes 
larger than one result in the model. These therefore should be regarded as upper 
limits. Finally the relatively simple energy models as well as the neglect of 
chemical production and loss terms are further simplifications which are very 
likely however masked by the uncertainties in the eddy diffusion coefficients and 
in the absolute content of atomic oxygen. 
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